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Abstract 

The connections between derivative operators and wavelets are well 
known. Here we generalize the concept by constructing multiresolu- 
tion approximations and wavelet basis functions that act like differ- 
ential Fourier multiplier operators. This construction follows from a 
stochastic model: signals are tempered distributions such that the ap- 
plication of a whitening (differential) operator results in a realization 
of a sparse white noise. Using wavelets constructed from these oper- 
ators, the sparsity of the white noise can be inherited by the wavelet 
coefficients. In this paper, we specify such wavelets in full generality 
and determine their properties in terms of the underlying operator. 



1 Introduction 

In the past few decades, a variety of wavelets that provide a complete and 
stable multiscale representation of L2(M'^) have been developed. The wavelet 
decomposition is very efficient from a computational point of view, due to 
the fast filtering algorithm. Underlying this approximation procedure is a 
differential operator, so our purpose in this paper is to construct wavelets 
that behave like a given differential Fourier multiplier operator L, which can 
be more general that a pure derivative. In our approach, the multiresolution 
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spaces are characterized by generalized B-splines associated with the opera- 
tor, and we show that, in a certain sense, the wavelet inherits properties of 
the operator. Importantly, the operator-like wavelet can be constructed di- 
rectly from the operator, bypassing the scaling function space. What makes 
the approach even more attractive is that, at each scale, the wavelet space 
is generated by the shifts of a single function. Our work provides a gen- 
eralization of some known constructions including: cardinal spline wavelets 
[5], elliptic wavelets fT5], polyharmonic spline wavelets [191 120], Wirtinger- 
Laplace operator-like wavelets [21j, and exponential-spline wavelets |13j . 

In applications, it has been observed that many signals are well repre- 
sented by a relatively small number of wavelet coefficients. Interestingly, 
the model that motivates our wavelet construction explains the origin of this 
sparsity. The context is that of sparse stochastic processes, which are defined 
by a stochastic differential equation driven by a (non-Gaussian) white noise. 
Explicitly, the model states that Ls = w where the signal s is a tempered 
distribution, L is a shift-invariant Fourier multiplier operator, and w is a 
sparse white noise [IB]. The wavelets we construct are designed to act like 
the operator L so that the wavelet coefficients are determined by a general- 
ized B-spline analysis of w. In particular, we define an interpolating spline 
0, corresponding to L*L, from which we derive the wavelets ip = L*(/>. Then 
the wavelet coefficients are formally computed by the L2 inner product 

{s,ij) = {s,L*(P) = {Ls,(P) = {w,^). 

Sparsity of w combined with localization of the B-splines results in sparse 
wavelet coefficients. This model is relevant in medical imaging applica- 
tions, where good performance has been observed in approximating func- 
tional magnetic resonance imaging and positron emission tomography data 
using operator-like wavelets that are tuned to the hemodynamic or pharma- 
cokinetic response of the system [12l [22] . 

Our construction falls under the general setting of pre-wavelets, which 
are comprehensively discussed by de Boor, DeVore, and Ron |7j. Two dis- 
tinguishing properties of our approach are its operator-based nature and the 
fact that it is non- stationary. Related constructions have been developed for 
wavelets based on radial basis functions [U EJ [17] . In fact [6] also takes an 
operator approach; however, the authors were focused on wavelets defined 
on arbitrarily spaced points. 

This paper is organized as follows. In Section [21 we formally define the 
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class of admissible operators and the lattices on which our wavelets are de- 
fined. In Section |3], we construct the non-stationary multiresolution analysis 
(MRA) that corresponds to a given operator L and derive approximation 
rates for functions lying in Sobolev-type spaces. Then, in Section HI we in- 
troduce the operator-like wavelets and study their properties; in particular, 
we derive conditions on L that guarantee that our choice of wavelet yields a 
stable basis at each scale. Under an additional constraint on L, we use this 
result to define Riesz bases of L2(M'^). In Section |5l we prove a decorrelation 
property for families of related wavelets. Finally, we conclude in Section [6] 
with some examples of the wavelets provided by our construction. 

2 Preliminaries 

The primary objects of study in this paper are differential Fourier multiplier 
operators and their derived wavelets. The operators we consider are shift- 
invariant operators L that act on the class of square integrable functions 
/ : M'^ — >■ C, which will be denoted by L2(M'^). The action of such a Fourier 
multiplier operator is defined by its symbol L in the Fourier domain, with 

L/= (l7)\ 

where the symbol L is a measurable function. Here / denotes the Fourier 
transform of / 

and (7^ denotes the inverse Fourier transform of g. Pointwise values of L 
will be required for some of our analysis, so we restrict the class of symbols 
by requiring continuity almost everywhere. Additionally, we would like to 
have a well-defined inverse of the symbol, so L should not be zero on a set of 
positive measure. To be precise, we define the class of admissible operators 
as follows. 

Definition 2.1. Let L be a Fourier multiplier operator, then L is admissible 
if its symbol L is a ratio of continuous functions f/g satisfying 

1. The set of zeros of fg has Lebesgue measure zero; 
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2. The zero sets of / and g are disjoint. 

Notice that each such operator defines a subspace of L2, consisting of func- 
tions whose derivatives are also square integrable, and our approximation 
results will focus on these spaces. 

Definition 2.2. An admissible operator defines a Sobolev-type subspace of 




Having defined the class of admissible operators, we must consider the 
lattices on which the multiresolution spaces will be defined. It is important 
to use lattices which are nested, so we consider those defined by an expansive 
integer matrix. Specifically, an integer matrix A, whose eigenvalues are all 
larger than 1 in absolute value, defines a sequence of lattices 

A^Z'^ = {A^k ■.keZ'^} 

indexed by an integer j. Here, we state several results about these lattices 
that will be used in this paper, most of which can be found in [11]. First, we 
know that A^U^ can be decomposed into a finite union of disjoint copies of 
^i+i^d. g^j^g |det(A)| vectors {e^jjfljj''^^'"^ such that 

IJ (A^e^ + A-'+^Z'^) = A^Z^ 
I 

and using this notation, our convention will be to set = 0. 

There are also several important properties that arise when using Fourier 
techniques on more general lattices. First, a lattice in the spatial domain 
corresponds to a lattice in the Fourier domain known as its dual lattice, and 
the dual lattice of A-'Z^ is given by 27r(A^)^-'Z"'. Also relevant is the notion 
of a fundamental domain, which for A-'Z"' is a measurable set i7j satisfying 

^ + A^fc) = 1 

for all X. Finally, just as in the dyadic case, we can define a discrete Fourier 
transform (DFT) corresponding to A in terms of the functions which are 
periodic with respect to the lattice. The matrix for this transform is given 
by 
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and the corresponding inverse DFT matrix is 

|det((A^)-i)r 

where k and I range over the cosets of 27c{A'^)^^'Z'^/'Z'^ and Z'^/AZ'', respec- 
tively. The inverse DFT matrix will play a role in deriving the Riesz basis 
condition on the wavelet spaces. We summarize the required properties in 
the following proposition. 

Proposition 2.3. A lattice generated by an expansive integer matrix A de- 
fines a discrete Fourier transform. This transform has a matrix representa- 
tion which is invertible and has nonzero entries. 



3 Mult iresolut ion Analysis 

The multiresolution framework for wavelet construction was presented by 
Mallat in the late 1980s p^. In the following years, the notion of pre- 
wavelets was developed, and a more general notion of multiresolution was 
adopted. We consider this more general setting in order to allow for a wider 
variety of admissible operators. 

Definition 3.1. A sequence {Vj}jQZ of closed linear subspaces of L2(M'^) 
forms a non-stationary multiresolution analysis if 

1. V,+i C V,-; 

2. \Jj^iVj is dense in L2(ffi'^) and f]j^^Vj is at most one-dimensional; 

3. f eVj if and only if /(■ - A^k) G Vj for all j G Z and k G Z'^, where 
A is an expansive integer matrix; 

4. For each j G Z, there is an element ipj G Vj such that the collection of 
translates {(fj{- — A^k) : k G Z'^} is a Riesz basis of Vj, i.e., there are 
constants < Aj < Bj < oo such that 



A, \\c\ 



^2 



< 



5^c[%,(--A^fc) 



< Bj \\c\\^^ 



L2 



5 



In order to produce non-stationary MRAs, we require additional prop- 
erties on the admissible operators. Together with a dilation matrix, the 
operator should admit generalized B-splines that satisfy decay and stability 
properties. 

As motivation for our definition, let us consider the one-dimensional ex- 
ample where L is defined by 

Lfit)^^{t)-am. 

In this Green's function for L is p{t) — e°*iJ(t), where H is the 

Heaviside function. In order to produce Riesz bases for the scaling matrix 
A = (2), we introduce the localization operator L^j defined by Ldj/ = 
/ — e^^"/(- — 2^). The generalized exponential-spline Lpj := L^jp will then be 
a compactly supported function whose shifts ipj{- — 2^) form a Riesz basis. 
Notice that, in the Fourier domain, is given by L~^{u)Ldj{u)) which is 



1 _ Q^^ia-iuj) 



lu — a 

In this form one can verify the equivalent Riesz basis condition 
Q<Aj<^ - 7r2-^'+^fe)f < Bj < oo. 



In fact, based on the symbol of L, we could have worked entirely in the 
Fourier domain to determine appropriate periodic functions L^ j. With this 
example in mind, we make the following definition. 

Definition 3.2. We say that an operator L and an integer matrix D are a 
spline-admissible pair of order r > d/2ii the following conditions are satisfied: 

1. L is an admissible Fourier multiplier operator; 

2. D = aR with R an orthogonal matrix and a > 1; 

3. For some constant cl > and all positive /i, we have 



1 + 



L{uj/h) 



-1/2 



Loo{^d.\Q.) 



where VL is the domain [— tt, tt]"'; 



6 



4. For every j G Z, there exists a periodic function L^j such that (pj{ijj) '■= 
Ld.j{uj)L{ij)~^ satisfies the Riesz basis condition 

0< Aj \i^j{uj + 27r(D^)--'fc)|^ < Bj < oo, 

kezd 

for some Aj and Bj in M. Here we require the periodic functions 
to be of the form ^kez^Pjl^^]^^^'^^'' some p G ii{Z'^). 

It can be verified that the last condition imphes that the function (pj is 
m L2(M'^). We call its inverse Fourier transform ipj := {<fj)^ a generalized B- 
spline for L. Also note that, in this definition, we have restricted the scaling 
matrices to be multiples of orthogonal matrices. The lattices generated by 
these matrices have some additional nice properties. For example, the lattices 
generated by D and are the same, and the lattices D-'Z'^ scale uniformly 
in j. Also, for such matrices, there are only finitely many possible lattices 
generated by powers of D; i.e., there always exists a positive integer n for 
which D" = a"I. 

Proposition 3.3. Given a spline- admissible pair L and D, the spaces 
form a non- stationary MRA. 

Proof. The fact that Vj+i C Vj follows from the definition of D and the 
Riesz basis conditions on the B-sphnes. The density property is a result of 
the admissibility of L, the Riesz basis condition, and the inclusion relation 
Vj+i ^ Vj, cf. [7\ Theorem 4.3]. In addition, the conditions imposed on L 
imply that the dimension of the intersection of the spaces in our MRA is at 
most one, cf. |71 Theorem 4.9]. Lastly, Properties 3 and 4 follow directly 
from the definitions. □ 

Condition 4 in the definition of spline-admissibility is the most challeng- 
ing to verify; however, there are many operators for which this condition is 
satisfied. For instance, it holds for any one-dimensional, constant-coefficient 
differential operator. In addition, it holds in higher dimensions for the Matern 
operators characterized by L{uj) = (1 + as they require no localiza- 

tion. In Section El we shall give a less trivial example and show how this 
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Riesz basis property can be verified. As a final point, note that if one is only 
interested in analyzing fine-scale spaces, Condition 4 need only be satisfied 
for j smaller than a fixed integer jo, but in this case, it is necessary to include 
the space Vj^ in the wavelet decomposition. 

We shall close this section by determining approximation rates for the 
spaces {Vj}j^z, in terms of L and the density of the lattice generated by D-'. 
In order to state this result, we must define the spline interpolants for the 

, 2 



operator L*L, whose symbol is 



The spline admissibility of this operator 



is the subject of the next proposition. 

Proposition 3.4. // L is spline admissible of order r > d/2, then L*L is 
spline admissible of order 2r > d. 

Proof. Let L be a spline admissible operator of order r > d/2. First notice 
that L*L is an admissible Fourier multiplier operator. Also, we can see that 
L*L satisfies Condition 3 of Definition 13.21 with r replaced by 2r. Therefore 
spline admissibility will follow if we can exhibit B-splines for L*L that satisfy 
the Riesz basis condition, where the integer dilation matrix is the same as 
for L. To this end, we claim that (pj := \(pjf defines the Fourier transform 
of a B-spline for L*L. 

An upper Riesz bound for ifj can be found by using the norm inequality 
between ii and £2- 

To verify the lower Riesz bound for (fj, we note that the decay condition 
of spline admissibility implies that for R > sufficiently large, we have 

Now for any M > sufficiently large, we can bound 

J2 |^j(w + 27r(D^)-^fc)|' 

above by using the decay estimate on L~^: 
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+ 27r(D^)-^fc)|' + C Ldj(u;) M'^-^' |det(D)|'^'"/'^ 



|fe|<M 

Therefore 



^ |^,-(a; + 27r(D^)-^fc)f > Aj - C L^,M 

\k\<M 



M'^-^'- \det(D] 



\2jr/d 



and 



|fe|<M 



|fe|<M 



>CM-'^ (^A, - C |Ld,j(u;)| |det(D)|'^'^/'^) 



Due to the fact that 2r > d, we can always choose M to make the right hand 
side positive, and this estabhshes a lower Riesz bound for ipj. □ 

The Riesz basis property of the B-splines for L*L imply that the L*L- 
spline interpolants 4>j{x), given by 



(oj) = |det(D) 



(1) 



are well-defined and also generate Riesz bases. Importantly, (pj G W2 does 
not depend on the specific choice of the localization operator, as we can see 
from 



0,H = |det(D)|^- 

= |det(D)|^' 





2 








L{u;) 



1 



1 + 
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These L*L-spline interpolants play a key role in our wavelet construction, 
which we describe in the next section; however, for our approximation result, 
we are more interested in the related functions 







-2 









-2 ' 



(2) 



which will also be needed for the decorrelation result in Theorem 15.41 

In order to bound the error of approximation from the spaces V^, we 
apply the techniques developed in [8j. In that paper, the authors derive 
a characterization of certain potential spaces in terms of approximation by 
closed, shift-invariant subspaces of L2(M'^). The same techniques can be 
applied in our situation, with a few modifications to account for smoothness 
being determined by different operator norms. 

The error in approximating a function / G L2(]R^) by a closed function 
space X will be denoted by 



£;(/,X) = min||/-s||j 

s&X 

and the approximation rate will be given in terms of the density of the lattice 
in W^. The lattice determined by D-' has density proportional to |det(D)p^'^, 
so we say that Vj provides approximation order f if for every / G VF2^(M'^) 



E(/,V,) <c|det(D; 



nf/d 



The rate of approximation is then determined by the following theorem. 

Theorem 3.5. For a spline- admissible pair L and D of order r > d/2, the 
multiresolution spaces Vj provide approximation order r < r if 



Idet(D^) 



-2jf/d 1 



1 + 



is bounded, independently of j , in Loo((D"^) ^Q), where Q = [— 7r,7r]'^. 

Proof. This result is a consequence of PI Theorem 4.3]. To show this let 
us introduce the the notation /j(-) = /(D-^-), which implies that fj = 
|det(D)|~-'/ o (D^)"-'. Notice that the spaces Vj are scaled copies of the 
integer shift-invariant spaces 
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We then can write the error of approximating a function / G 14^2^ (M'') from 
Vj in terms of approximation by V- as 



E{f,V,) = Idem f'E{f,,V^^) 

= (27r)-^/2|det(D)r/^E(g,V7-), 



where Vj is composed of the Fourier transforms of functions in Vj . Separat- 
ing this last term, we have 



E{f, V,) < {2n)-'/' |det(D)p/^ iE{f,xn, K^) + (1 - Xn)f: 



li/2 



(3) 



where Xfi is the characteristic function of the set Q. We are now left with 
bounding both terms on the right-hand side of ([3]). First, we have 



[^-Xn)fj 



duj 



|det(D) 



-2j 



1 + 



L((D^)-^u;; 



1 + 



L{{-DT)-Juj) 



and since L is spline- admissible of order r 



;i - Xn)f, ^ < |det(D; 



i2jr/d-j II f II 2 



SO that 



|det(D)|^-/^ il-Xn)f, 



< IdetfD' 



nr/d 



(4) 
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In order to bound the remaining term, we will need a formula for the projec- 
tion of fjXn onto Vj . Notice that 



1 -mj((D^) -■'■•) = 1 



Efc«.|^.°(D^)-^'(--27rfc)r 

(fj o 



so we can apply ^ Theorem 2.20] to get 



2 ' 



;i-m,((D^)-^-.)) 



= |det(D)r'^' 
Now, changing variables gives 



/o(D^)-^- (l-m,((D^)-^-.)). 







2 




1 


f 


(1 + 


L 











2 1 — rrij 



1 + 



< IdetfD) 



-i II /II 2 



1 - 


m 


j 






2 


1 + 


L 





Applying our assumption on (1 — mj), we have 

\demf'E{f,xn,V^n<c\demr/'lU\\w- 
Substituting the estimates and ([5]) into (E]) yields the result. 



(5) 
□ 



4 Operator-Like Wavelets and Riesz Bases 

Using the non-stationary MRA defined in the previous section, we define the 
scale of wavelet spaces Wj by the relationship 
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i.e., VTj+i is the orthogonal complement of V^+i in Vj. Our goal in this section 
is to define Riesz bases for these spaces and for L2(]R'^). To begin, let us define 
the functions 

which we claim generate Riesz bases for the wavelet spaces, under mild con- 
ditions on the operator L. First, note that i^j+i is indeed in V}, because its 
Fourier transform ■0^+1 is a periodic multiple of (pj, and thus 

V^,+i(u;) = |det(D)|^' L^A^)* ^ 

A direct implication of our wavelet construction is the following property. 

Property 4.1. The wavelet function ■^j+i behaves like a multiscale version 
of the underlying operator L in the sense that, for any / G VF^, we have 
/ * ^j+i = L*{/ * <Pj}- Hence, in the case where (j)j is a lowpass filter, 
{L*{/ * 4>j}}j£Z corresponds to a multiscale representation of L*{/}. 

The next few results focus on showing that the D-'Z'^/D^+^Z'^ shifts of 
are orthogonal to Vj+i and generate a Riesz basis of Wj+i. 

Proposition 4.2. The wavelets — 'D-'k)}i^^j,d\pzd are orthogonal to 

the space V^+i- 

Proof. Note that it suffices to show (</7j+i, — D-'fc)) = for every 

k e Z'^\DZ'^. By definition 

(^,+i,^,+i(--D^fc))= / ^,+i(u;)e^'^-°^'=L(u;)0,(a;)du; 

= / Ld,,+i(u;)e^-^^*'J,(u;)du;. 

Now let Q be a fundamental domain for the lattice generated by 27r(D^)~-' 
so that we can write 
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= 0. 

The last equality holds because i^d,i+i is periodic on a finer lattice. □ 

In order to prove that shifts of ipj form a Riesz basis of the wavelet space, 
we need to introduce some notation, which will help formulate the problem as 
a shift-invariant one. Since we view the wavelet space as being shift-invariant 
on D-'"*"^Z'^, a representative wavelet is chosen for each of the |det(D)| — 1 
cosets. Specifically, recall that there are vectors {e;}|jfQ*'^^'~^ in Z*^ such that 

[j [B^ei + D^+^Z'^) = D^Z^ 
so for 1 < / < |det(D)| — 1, we define the wavelets 

In the following, necessary and sufficient conditions on the operator L will 
be given which guarantee that \1/ = generates a Riesz basis 

of Wj+i. The technique used is called fiberization and it can be applied to 
characterize finitely generated shift-invariant spaces [16]. In this setting, a 
collection of functions defines a Gramian matrix, and the property of being 
a Riesz basis is equivalent to the Gramian having bounded eigenvalues. In 
our situation, the Gramian for \l/ is 




-iDJ(ej.-e;)-{a)+a) 



k,l 



Let us denote the largest and smallest eigenvalues of G^r by A (a;) and 
A (a;), respectively. Then the collection \1/ is a Riesz basis if and only if A 
and 1/A are essentially bounded (cf. [12] Theorem 2.3.6). To simplify this 
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matrix without changing the eigenvalues, we apply the similarity transfor- 
mation T{u})^^G^l,{u:)T{(JJ), where T is the diagonal matrix with diagonal 
entry e~*°^'^' '^ in row /. This transformation multiplies column / of Gq; by 
g-jD^ero; ^ q£ giDJefe W gjj^ce the eigenvalues are unchanged, 

let us call this new matrix as well. After applying the transformation, 
we have 



^ae27r(DT)-3-iZd 

Using the fact that |J^ (e^ + Ti^l/^ = and the notation 
c(m;u;) := J] + 27r(D^)-^'-^e^ + /3) 

/3627r(DT)-3Zd 

we can write 

Vet(D)|-l 

6-^(0;)= I J] c(m; u;)e-''^'(^'=-^'>(°^)"'^" 



m=0 



Since is Hermitian, its maximum and minimum eigenvalues can be de- 
termined by analyzing the quadratic form 



a 



|det(D)|-l 

'G^OL = Y ocloci ^ 

k,l=l 
|det(D)l-l 



det(D)|-l 

c(m; u))e~ 

m=0 

|det(D)|-l 



■2m{ek-ei)-(D'^)-^e„ 



= J2 c(m;a;) {ake^-'-'^<''^r'emy^^^2nier{n-)-^er, 



m=0 



|det(D)|-l 

J2 c(m;uj) 



m=0 



k,l=l 
|det(D)|-l 



J2 a^e^--^-^^' 



k=l 



over all unit vectors a e Cl'^^**^^)! ^. Let us additionally define H to be the 
inverse DFT matrix for the lattice generated by D, i.e.. 
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V / k,m. 

with k and m ranging from to |det(D)| — 1, and let Hq be the submatrix 
obtained by removing row from H. Note that each sum in the quadratic 
form above is given by the inner product of a with a column of Hq. This 
identification allows us to prove the following. 

Lemma 4.3. For a fixed u:, the quadratic form cx.*Gq,a is non-zero if and 
only if at most one of the coefficients c{m; cv) is zero. 

Proof. First, suppose that two distinct coefficients, c(mo; a;) and c(mi;a;), 
are zero. Then we can clearly find an ck that is orthogonal to the remaining 
|det(D)| — 2 columns of Hq. Hence the quadratic form will be for such an 
ct. 

Next, suppose that all coefficients except possibly c(mo; uj) are non-zero, 
and denote the matrix obtained from Hq by removing column mo by Ho,mo- 
Then, to prove positivity of the quadratic form, we show that the columns of 
Ho,mo span C''^'^'^^^'^^. This is equivalent to showing that the determinant of 
Ho,mo is non-zero, which follows from the fact that has non-zero entries 
(since Ho,mo is a first minor of H). □ 

Lemma 4.4. The collection is a Riesz basis if and only if no two of the 
coefficients c{m; u:) are zero for the same uj. 

Proof. This follows from the definition of and Lemma 14.31 □ 

Suppose that the conditions of Theorem 14.51 are satisfied so that is 
a Riesz basis. Then the orthogonality between the shifts of \E' and V^+i, 
together with the Riesz basis condition on v'j+i, implies '^'■.={ipjj^i\ U \l/ 
generates a Riesz basis. The fact that spans Vj follows from a comparison 
with the Riesz basis generated by the D-'+^Z'^ translates of ipj{- — D-'e;), cf. 
[1] and [TJ Theorem 2.26]. As a consequence of these results, ^ would be a 
stable basis of Wj+i. 

In order to interpret this result in terms of the operator L, we note that 
'tpj^i^oj) vanishes at p + (D'^)~-^27rfc for every p & such that L{p) = 0. 
In other words, each zero of L generates a periodic sequence of zeros in the 
spectrum of the generating wavelet. As long as these periodizations of zeros 
are essentially disjoint, the conditions of Theorem 14.51 will be satisfied. The 
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zero set of the symbol is denoted hj N = {p E R'^|L(p) = 0}, and for each 
scale j, we define 

ATf = p + (D^)-(^'+i)27rei + (D^)-^27rZ^ / = 0, . . . , |det(D)| - 1, 
where p G M . 

Theorem 4.5. Let j E 7j be an arbitrary scale. Then the family of functions 
^ = ^ generates a Riesz basis of Wj+i if and only if the sets 

Mf^ satisfy 

for each 1 < / < |det(D)| - 1. 

Proof. This follows from the periodicity of the zeros of ^j+i and the condition 
on the coefficients from Lemma I4.4[ □ 

Our wavelet construction is intended to be general so that we may account 
for a large collection of operators. As a consequence of this generality, we can 
not conclude that our wavelet construction always produces a Riesz basis of 
L2(M'^). Here, we shall impose additional conditions to ensure a Riesz basis 
is produced. In order to preserve generality, we shall focus on the fine scale 
wavelet spaces and include a B-sphne space V^^ in our Riesz basis. 

Theorem 4.6. Let L be a spline admissible operator of order r, and suppose 
that there exist uo > and constants Ci,C2 > such that the symbol L 
satisfies 



Ciiu^r < 



Liu;] 



< Co \uj\ 



for \uj\ > Wq. Then there is an integer jo such that the collection 

- m^m- U {|det(D)|(^'+^)/(''^)-^' - /3)}^,D.z^\D.+iz^ (7) 

j<jo 

forms a Riesz basis of L2(M.'^). 
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Proof. We shall verify the result for the dilation matrix D = 21, as the other 
cases are similar. Let jo be an integer for which cjq < 27r2^-'"~^. 

Considering Lemma l47i| we can see that the Riesz bounds for the wavelet 
spaces depend on estimating the functions c(m; u), the main term being the 
sum 



Y^^^^a + 27r2-i-ie„ + f3 + 2n2-ik) 
The value of this sum depends on the position of the lattice 



-2 



(8) 



Xj{m,u:] 



u + 27r2-^'-^ 



27r2-^Z'^ 



with respect to the origin. To prove bounds of ([HD, let us introduce the 
notation hj for the fill distance and qj for the separation radius of the lattices 
Xj{m,u})] i.e, 



hj := sup inf \x — Xk\ = 27r2 ^ ^Vd 



and 



inf \xk — Xk'\ = 27x2 ^ ^. 

k^k' 



Since each lattice Xj{m, uj) is a translation of Xj{0, 0), the quantities hj and 
Qj are independent of m and 

We shall now bound (|8]) by considering two cases: 

1. dist{0,Xj{m,u})) > qj/2- 

2. dist(0,Xj(m,Lj)) < qj/2. 

In the first case, all points of the lattice Xj(m, cj) lie outside of the ball 
of radius ujq centered at the origin. Therefore, ([H]) can be reduced to 



5^ \L{xk, 



(9) 
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Applying Proposition lA.ll and Proposition \A.2\ we can find upper and lower 
bounds for this sum which are proportional to 2^^~^. Importantly, the pro- 
portionality constants are independent of j. 

For the second case, we must be more careful as one of the lattice points 
will lie close to the origin. For any fixed u, there will be at most one m 
for which the lattice Xj(m, oj) satisfies this second condition. Therefore, a 
sufficient lower bound for the lattice sum will be 0; however, the upper bound 
must match the one derived above. Let us consider the cases 

2a. L takes the value on the lattice Xj{m, oj) 



2b. 



L- 



< oo on the lattice Xj{m^ oj) 



If the first condition is satisfied, then ([H]) will be 0. However, if the second 
condition holds, then we can again reduce the sum to ([9]). Clearly, is a 
lower bound for this expression. For the upper bound, we apply Proposition 
lA.ll as in the previous case. 

To finish the proof, we note that Lemma 14.31 implies that the wavelets 

form a Riesz basis at each level j < jo- Furthermore, the bounds on c(m, uj) 
imply that the Riesz bounds are all equal. Therefore, the collection ([7]) is a 
Riesz basis of La (M'^). □ 



5 Decorrelation of Coefficients 

As was stated in the introduction, the primary reason for our construction 
is to promote a sparse wavelet representation. Our model is based on the 
assumption that the wavelet coefficients of a signal s are computed by the 
L2 inner product 

(z^,(--D^fc),s>. 

Here, we should point out that, unless the wavelets form an orthogonal basis, 
reconstruction will be defined in terms of a dual basis. However, as our focus 
in this paper is the sparsity of the coefficients, we shall be content to work 
with the analysis component of the approximation and leave the synthesis 
component for future study. 
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Now, considering our stochastic model, it is important to use wavelets 
that (nearly) decorrelate the signal within each scale, and one way to ac- 
complish this goal is by modifying the underlying operator. Hence, given a 
spline-admissible pair (L,D), we define a new spline-admissible pair (L„,D) 
by L„ := L", and we will see that as n increases, the wavelet coefficients be- 
come decorrelated. This result follows from the fact that the (L„)*L„-spline 
interpolants (appropriately scaled) converge to a sinc-type function, and it 
is motivated by the work of Aldroubi and Unser, which shows that a broad 
family of spline-like interpolators converge to the ideal sine interpolator [2J. 
To state this result explicitly, we denote the B-splines for L„ by (pn,j = 
Therefore the (L„)*L.„-spline interpolants are given by 



(a;) = |det(D 

= |det(D)p 
and we analogously define 



E.ex.|^n,,(a; + 27r(D^)-^-fc)|^ 



E.ez^l^.(^ + 2vr(D^)-^-fc)| 



2n ' 



Decorrelation will depend on showing that the functions m„ j converge to 1 
almost everywhere on some fundamental domain VLj of 27r(D^)~-'Z'^ and to 
almost everywhere outside of Vtj. Our proof relies on techniques used in 
the analysis of radial basis functions. For example, a similar method was 
used by Baxter to prove convergence of the Lagrange functions associated 
with multiquadric functions [3", Chapter 7]. The idea is to define disjoint sets 
covering R'', each with a single point in any given fundamental domain, and 
analyze the convergence of m„j- on these sets. 

Definition 5.1. For each j G Z and x G Vtj we define the set 

Ej^^ := {x + 27r(D^)~^fc : k G Z'^}. 

By the Riesz basis condition, each Ej^x has a finite number of elements y 
with mij{y) of maximal size, and we define Fj to be the set of x G Qj such 
that there is not a unique y G Ej x where rrii j attains a maximum. 
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Lemma 5.2. Let x e i}j\Fj, then for y e Ej ^ we have rrinjiy) — )■ «/ and 
only if mi j{y) is not of maximal size over Ej ^- Furthermore, if mij{y) is 
of maximal size, then runj {y) 1 . 

Proof. Fix X G ^j\Fj and y e Ej^^- Notice that the periodicity of the 
denominator of rrii j impUes that mij{y) is maximal iff \(fij{y) \ is maximal 
Let us first suppose mi_j{y) is not maximal. If \(pi^j{y) \ = 0, the result is 
obvious. Otherwise, there is some ko G Z'^ and 6 < 1 such that 

Therefore 



and the result follows. 

Next, suppose mij{y) is of maximal size. Since (pij has no periodic zeros, 
\^i,j{y) \ 0- Therefore 



TT^n,j{y) = 



1 



BnAy) 



with 



BnAy) = E 



My) 



Since \^Pi,j{y)\ is of maximal size, all terms of the sum except one are less 
than 1. In particular, B^j will converge to 1 as n increases. □ 

Lemma 5.3. Given j e Z, if the Lebesgue measure of Fj is 0, then 



converges to 1 in Li{Qj) as n — >■ oo. 
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Proof. Notice that the sum is bounded by 1 . Therefore Lemma 15.21 imphes 
that the sum converges to 1 on the complement of Fj. Hence, we can apply 
the dominated convergence theorem to obtain the result. □ 

With this theorem, we can show how the wavelets corresponding to L„ 
decorrelate within scale as n becomes large. The way we characterize decor- 
relation is in terms of the semi-inner products 



{f,9)r 



fg* 



which is a true inner product for the wavelets 

Theorem 5.4. Suppose the Lebesgue measure ofUj^zF'j, as provided in Def- 
inition [377]. is 0. Then as n increases, the wavelet coefficients decorrelate in 
the following sense. For any j ^ Z,k & Z'^\{0} we have 

(V'nJ+l, V'nj + l(- - -> 

as —7- oo. 

Proof. The inner product can be expressed as an integral using the definition 



-2 



da; 



|det(D)|'^' / mn,j{uj)^e-'^"''^duj, 



and we can periodize the integrand to get 



m„j(LJ) e 



da; 
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The last expression converges to by the Lebesgue dominated convergence 
theorem. □ 



Let us now show how this resuh imphes decorrelation of the wavelet 
coefficients. To begin, we introduce the notation 

Cn,j+l,k ■= {i^n,j{- - D-'fc), s) 

for the wavelet coefficients of the random signal s. Note that our stochastic 
model, Ls = w, implies that the coefficients are random variables. Hence, for 
distinct k and k', the covariance between Cnj+i^k and Cnj+i^k' is determined 
by the expected value of their product: 

E{cnj+i,fcC„j+i,fe'} = E{(0„j(- - I)^k),w) (0„j(- - D^fc'),w)} . 

Now, as long as the white noise w has zero mean and second-order moments, 
the covariance satisfies 

E{Cn,j+l,kCn,j+l,k'} = {(t>n,j{- " D'' fc) , j (■ - D^fc')) 

where the convergence follows from Theorem 15.41 Therefore, when the stan- 
dard deviations of c„,j+i,fc and Cnj+i,k' are bounded below, the correlation 
between them will also converge to zero. 

6 Examples 

The wavelet construction that we have presented in this paper is quite general 
and accommodates many operators. In this section, we consider operators 
that are defined as polynomials in the Laplacian. In particular, we shall 
show how spline admissibility may be verified for such operators. Also, note 
that each of these operators satisfies the growth condition of Theorem 14.61 
Therefore, the corresponding wavelet spaces may be used to construct Riesz 
bases of L2(M'^). 
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6.1 Matern and Laplace Operators 



First, consider the two-dimensional Matern operator that is defined by the 
symbol Ly{uj) = (|cl»|^ + 1)''/^, with the parameter u > 1. As L^{(jj)~^ satisfies 
the Riesz basis condition, no localization operator is needed. Therefore, the 
operator L^, is spline-admissible of order u. Also note that the fractional 
Laplacian L = \u:f and iterated Laplacian L = \ujf operators are also 
spline-admissible. 

Each of these wavelets can be considered with any admissible subsampling 
matrix D. In two dimensions, the choice of D = 21 corresponds to the 
rectangular grid, while 

characterizes the quincunx subsampling scheme. 



6.2 Helmholtz Operator 

Operators whose symbols have more general zero sets can also be used. As 
an example, consider the polynomial p{x) = 1/4 + x. The operator p{A) is 
known as the Helmholtz operator, and its symbol is L^cu) = 1/4 — \ujf. The 
wavelets corresponding to this operator could potentially have applications 
in optics, as the Helmholtz equation. 



An + Am = /, 

is a reduced form of the wave equation [9l Chapter 5]. In what follows, we 
show that this is a spline-admissible operator for the scaling matrix D = 21 
on R^. However, since the wavelets ipj = L*0j will not form a Riesz basis for 
the coarse-scale wavelet spaces, we only consider j < —1. 

The primary issue with proving spline-admissibility is defining discrete 
operators Ldj corresponding to L. It is known that sufficiently smooth func- 
tions have absolutely convergent Fourier series [TOj Theorem 3.2.9], so we 
must construct such a function with the same zero set as L. Since L is zero 
on a circle and radial, we can choose L^ j to be the periodization of a radial 
function. More precisely, we define fj to be a function that is periodic with 
respect to the lattice dual to 2-'Z^ and, on [— 27r,27r]^, satisfies: 
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1. The function fj{oj) has continuous derivatives up to order 3 and is 
constant for |a;| > 1; 

2. There are constants Ci and C2 such that 



where 02,fc(M) is a Wendland function with 2k continuous derivatives [23^ Sec- 
tion 9.4]. We can therefore apply our wavelet construction, which produces 
a Riesz basis. 

A Discrete Sums 

Let X = {xk} be a countable collection of points in M'^, and define 



Also, let B{x,r) denote the ball of radius r centered at x. Proving Riesz 
bounds for the wavelet spaces relies on the following propositions concerning 
sums of function values over discrete sets. 

Proposition A.l. If hx < 00 and r > d/2, then there exists a constant 
C > ( depending only on r and d, not hx ) such that 



< Ci < fj{uj){l/4: - \u)\^) < C2 < cx) 



for |a;| < 1. 



These conditions are met, for example, by defining fj by 





\xk\ > Ch 



-2r 

'X 



Xk\>2hx 



Proof. Since x^ > 2hx, we have 



Xk - hx 



> 2 



-1 
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which imphes 



Then 



Xk-h 



X- 



-2r 



E 



|-2r ^ I |-2r 

Xk\ > 2^ \Xk\ 



Xk>2hx 
-2r 



>2- Yl 

■k 
2 



2r 



Xk - hx 



Xk 



Xk 



Vol {B{xk,hx)) 
Vo\{B{xk,hx)) 



> 



E 



Vol(B(0,Ax))^^t;i^yB(,.,»,) 



/•oo 



I3hx 



□ 



Proposition A. 2. // \xk\ > ?x/4 and r > d/2, then there exists a constant 
C > ( depending only on r and d, not qx ) such that 



J2 k^r"^ < cq-, 



2r 



Proof. Using the fact that |£Cfe| > gx/4, we can write 

qx Xk 



Xk + 



8 \xk\ 



< 2 \xk\ , 



which imphes 



Xk + 



qx Xk 



8 \xk\ 



-2r 



We now have 
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< 



Xk 
22r 



qx Xk 



Xk 



E 



Vol{B{xk,qx/8)) 
Yo\{B{xk,qx/8)) 



< 



Vol(i?(0,gx/8))^^(.„,,/8) 

22r 



\x\-^''dx 



Vol(S(0,gx/8))^|>,,/8 

POD 

Jqx/S 

< Cqx''- 



x\ dx 



□ 
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